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We study current induced magnetization dynamics in a long thin ferromagnetic wire with 
Dzyaloshinskii-Moriya interaction (DMI) . We find a spiral domain wall configuration of the magne- 
tization and obtain an analytical expression for the width of the domain wall as a function of the 
interaction strengths. Our findings show that above a certain value of DMI a domain wall configu- 
ration cannot exist in the wire. Below this value we determine the domain wall dynamics for small 
currents, and calculate the drift velocity of the domain wall along the wire. We show that the DMI 
suppresses the minimum value of current required to move the domain wall. Depending on its sign, 
the DMI increases or decreases the domain wall drift velocity. 

PACS numbers: 75.78.Fg; 75.60.Ch; 71.70.Ej 



Introduction. A number of recent experiments, per- 
formed in various metallic magnets, have shown the spiral 
structure of magnetization due to Dzyaloshinskii-Moriya 
interaction (DMI) [1, 2, 3, 4, 5, 6, 7, 8]. In particular, 
the B20 structure of ferromagncts, such as MnSi, which 
lacks strict space-inversion symmetry, leads to a long- 
wavelength helical twist in the magnetization [4, 5, 6]. 
Furthermore, the direct space-time observation of the spi- 
ral structure by Lorentz microscopy became possible for 
DMI- induced helimagnets [7, 8]. Using spin-polarized 
scanning tunneling microscopy it has been shown that 
the magnetic order of 1 monolayer Mn on W(001) is a 
left-handed spiral [2] and that the magnetic structure of 
the Fe double layer on W(110) is a right-rotating spiral 
[1]. All these spiral states are the consequence of DMI. 

A spin-polarized current flowing through such spiral 
magnetic structures would exert a spin-torque which 
could be used for manipulations of the magnetization 
with potential applications. For example, in magnetic 
memory devices [9, 10] the key issue is to manipulate the 
domain wall (DW) configurations by means of magnetic 
fields and/or spin-polarized current. Therefore, current- 
induced dynamics of spiral magnets is an important sub- 
ject of technological relevance. 

One of the most important factors which effects the 
DW motion is pinning. The DW pinning can have "ex- 
trinsic" and "intrinsic" nature. The extrinsic pinning 
is due to surface roughness and other irregularities of 
the wires which brake translational invariancc. On the 
other hand, the intrinsic pinning is present even in ideally 
smooth (translation invariant) nanowires. It depends on 
the wire geometry and material parameters which can 
be described by anisotropics. Although extrinsic pinning 
can be significantly reduced in the near future with the 
help of more sophisticated wire fabrication techniques, 
the intrinsic pinning is always present. Therefore, in this 
Letter we concentrate on the more important case of DW 
dynamics with the intrinsic pinning. 

We determine the effect of a polarized current on the 



magnetization configuration in the ferromagnetic wire 
with both strong easy-axis anisotropy along its axis and 
weak anisotropy in the plane transverse to the wire. The 
DMI, which arises from spin-orbit scattering of electrons 
in non-centrosymmctric magnetic materials is typically 
irrelevant in bulk metals as their crystals are inversion- 
symmetric. However, in low-dimensional systems (such 
as atomic layers and wires), which lack structural inver- 
sion symmetry, the DMI in the presence of softened fer- 
romagnetic exchange coupling leads to the formation of 
the spiral spin structures. 

The main goal of this Letter is to study the influence 
of DMI on the magnetization dynamics in ferromagnets. 
We obtain the expression for the DW width as a function 
of the DMI constant , uniaxial anisotropy along the wire, 
and exchange interaction constant. Wc find that there 
is a critical value of the DMI above which a DW config- 
uration cannot exist in the wire. This result can have 
an important implication for the future experiments by 
setting a limit on the devices with DMI which use DWs 
for information manipulation. Below this critical value 
of DMI the DW can propagate along the wire and ro- 
tate around its axis. Any angle is equally favorable for 
the DW if there is no anisotropy in the transverse plane. 
Generally speaking in most wires there exists such an 
anisotropy due to the asymmetry of the wire cross sec- 
tion. We show that it leads to a chosen direction of the 
magnetization in the center of the DW, so that the wall 
cannot rotate freely anymore. Therefore, if a polarized 
current is passing through such a wire, the DW will move 
only if the current is larger than a certain critical value. 
This value corresponds to the minimal torque needed to 
be pumped into the system to rotate the spins of the DW 
around the wire's axis. 

We investigate the dynamics of the DW in the small 
transverse anisotropy regime. In particular, we find the 
drift (average) velocity of the DW in the wire with DMI. 
Our findings also show that DMI decreases the critical 
value of current required to move a DW. To obtain all 
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these results for the DW dynamics, we use a universal 
method for finding zero mode dynamics of spin textures. 
This method is described in detail in the supplementary 
material [11]. 

Model. We employ a simple theoretical model of a fer- 
romagnet with DMI and anisotropics which highlights a 
new kind of behavior of DW structures. We consider a 
Hamiltonian for a fcrromagnct which has two terms de- 
scribing the exchange and Dzyaloshinskii-Moriya inter- 
actions [12, 13]. Without anisotropics in the continuous 
limit this Hamiltonian takes the form. 
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A)M • (V x M) 



(1) 



Here M is a magnetization vector, Jo > is exchange 
interaction constant, and Do is the DMI constant. Wc 
study a ferromagnetic wire which is modeled as a one- 
dimensional (ID) classical spin chain [21], where the wire 
is along the z-axis, see Fig. 1. For the thin long wire with 
uniaxial anisotropy Hamiltonian (1) modifies to 
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FIG. 1: (color online) Sketch of the wire with magnetization 
profile for a DW. A and K denote the anisotropies along and 
transverse to the wire, respectively. The upper inset shows 
the dependence of j c on the twist FA, Eq. (18); the lower 
inset shows drift velocity Vd and variance ((SV) 2 ) T (in arb. 
units) vs. current j, see Eqs. (20) and (21). 



Here e 2 is the unit vector in z direction, d = d/dz, 
and we introduced normalized magnetization vector S = 
M/M, so that S 2 = 1, D = D /(AM 2 ), and J = 
Jo / {AM 2 ), where A is the cross-sectional area of the 
wire. The last term in Eq. (2) is due to uniaxial 
anisotropy (with the anisotropy constant A = Ao / '(AM 2 )) 
which shows that the system favors the magnetization 
along the wire. 

To study the magnetization dynamics we employ the 
generalized Landau-Lifshitz- Gilbert equation [14, 15] for 
ID wire with current j along the wire: 



S = S x H e - jdS + 0jS x as + aS x S. 



(3) 



where H e = 5H/5S, S = dS/dt, a = a^/M 2 and ao is 
the Gilbert damping constant, f3 = fto/M 2 and @o is the 
constant of nonadiabatic current term, time is measured 
in the units of the gyromagnetic ratio 70 = g |e| /(2mc), 
and the current j is measured in units of a 3 / (2eM7o) 
where a is the lattice constant. Generally speaking one 
also has to specify the boundary conditions for Eq. (3). 

A general solution of one-dimensional LLG equation 
(3) can always be presented in the form 



as = T(z, t)e z x S + A(z, i)Sx [e z x S] , 



(4) 



where T and A are in general two independent functions 
of z and t; it also follows that dS z = A(l — S 2 ). 

Zero current. First we consider the simplest case of 
zero current (j = 0) and look for a time-independent 
magnetization configuration. This means that we need 
to minimize Hamiltonian (2) which can be written up to 
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The spin configuration depends on the sign of X — D 2 /2J. 

For 2 JA < D 2 the minimum of the second term is at 
S z = 0. The first term is minimized by the condition 
as = -je 2 x S, so that the solution is a spiral, 

S = (cos(r2 + 0o),sin(rz + </>o),O) T , T = D/J. (6) 

The ground state is thus unique and there is no DW 
configuration. Therefore, for the wires with weak enough 
uniaxial anisotropy and/or exchange constant the spiral 
magnetization state can prevent the formation of DWs. 

For 2 JA > D 2 the minimum of the second term is at 
S z = ±1. This also minimizes the first term in Eq. (5). 
Thus, S z = ±1 are the two solutions, and a DW can exist 
in the wire as a transition from one solution to another. 

Below we study the statics and dynamics of such a 
DW in the wire, and therefore we concentrate on the 
case 2JA > D 2 . Then the boundary conditions for Eq. 
(3) are S z — ► ±1 at z — ► ±00. To find the static configu- 
ration of the DW we consider the solution in the form (4). 
Substituting it into Hamiltonian (5), we find 
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The minimization of the first term sets 
T = D/J, 



(7) 



(8) 
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cf. Eq. (6). Using paramctrization S z = tanh/(z), we 
obtain 

(a/) 2 + A~ 2 



H = 



J 



dz- 



cosh / 



(9) 



where A 2 , = <J J/2A is the DW width in the absence of 
DMI. The straightforward minimization of Eq. (9) gives 
/ = z/A or A = 1/A, and in components the solution 
takes the form 



= cos(r(z - go) + 4>) 

cosh((z — zo)/A) 
= sin(r(z - z ) + <f>) 
y cosh((z- z )/A) ' 
S z = tanh(z — zq)/A), 



(10a) 

(10b) 
(10c) 



where the angle <fi is the tilt of the DW, and zq is its 
position (both arbitrary). We see that 2tt/T is the pitch 
of the spiral, A is the width of the DW and T A is thus the 
twist of the DW. Both T and A have the same functional 
dependencies in terms of Jo, -Do, and Ao as in terms of 
J, D, and A. According to its definition in Eq. (9), A 
becomes infinite in the limit 2JA = D 2 and DW cannot 
be sustained in the wire. 

The energy of the DW is E = 2J/A = 2^2J\- D 2 . 
It vanishes when D 2 approaches 2JA. 

The z component of the magnetization (10) is the same 
as that of a standard (without DMI) DW of width A. 
The direction of the twist of the DW depends on the 
sign of DMI and can be either clock or counterclockwise. 

Parameters zq and (f> in Eq. (10) correspond to two 
zero-modes of the system. These modes are the most 
relevant if the system is perturbed. The time-dependent 
solution then can be represented in the form of a moving 
and rotating DW plus a small correction to its shape. 
The requirement that the correction remains small dur- 
ing the motion leads to the equations for the velocity 
and angular velocity of the DW. A detailed derivation of 
these equations is presented in the supplementary mate- 
rial [11]. Below we present the results and discuss their 
implications. 

Small currents. First we find the magnetization dy- 
namics in the wire for small applied currents. We denote 
the solution (10) for the DW without a current, as So(z). 
When the current is applied we expect the DW to move 
and rotate. The full dynamics is described by the equa- 
tion 

S = S x H e + h, H e = SH/SS, (11) 

where the correction h for small currents is h = hy , 

hj = -jdS + ((3j-az )S xdS + a(j)S xe z xSo. (12) 

This correction gives the following results for the DW 
velocity and angular velocity: 

• l + aP+ja- g)TA . • (a - /?)A . 



A few conclusions can be made from these equations. 

i. ) The direction of the DW rotation depends only on 
the relative strength of the two dissipative terms in the 
LLG Eq. (3). Remarkably, the sign of the DMI correction 
to the DW velocity depends on weather the DW rotates 
in the same direction as the twist of the DW. 

ii. ) For p = there is an integral of motion a(r 2 + 
A~ 2 )z - (1/A + aT)(j) = const. If we take 1/A = 
which corresponds to a perfect spiral state (DW width is 
infinite), this invariant just describes the rotation of the 
spiral while it moves. 

iii. ) At very large twists TA, zo = 4>TJ/2X, indepen- 
dently of both a and (3. 

iv. ) The DW rotation and its velocity diverge when 
D 2 approaches 2A J. This nonphysical result is the con- 
sequence of the fact that our derivation of Eq. (13) 
neglects all modes except the zero ones. However, in the 
limit of D 2 — > 2AJ the breathing mode (the mode that 
corresponds to the change of the DW width and pitch) 
softens and its dynamics cannot be neglected [16]. 

In line with the general result [17], in the special case 
of a = (3, the DW docs not rotate and just moves with 
the velocity which depends on current only. 

Small anisotropy in the transverse plane. In order to 
account for the anisotropy in the transverse plane we 
introduce a correction to Hamiltonian (2) in the form 



dzKS 2 (z) 



(14) 



where the anisotropy constant K > is very small. 

The presence of this anisotropy fixes the tilt angle <j> of 
the solution (10). To show it we calculate the correction 
to the energy to the first order in K by substituting Eq. 
(10) into Eq. (14). Assuming the wall to be at the origin 
(z = 0), we obtain 5 X E = KA - B ?* h ^) c os{24>). This 
correction has a minimum at </> = 0, ir. When DMI is ab- 
sent (r = 0) this correction reduces to K A[l — 2cos(2</>)]. 

Dynamics and transverse anisotropy. Now we find how 
the small anisotropy in the transverse plane affects the 
magnetization dynamics. The correction h defined in 
Eq. (11) takes the form h = hj + h xy , where hj is given 
by Eq. (12) and 

h xy = S x = 2KS y S x e y . (15) 

This perturbation leads to the following equations for 
the position of the DW and the tilt angle: 

(3 (a- /3)(l + aTA) . . ,„ , N1 , „. 

Zo = * + 7TZ 2^ " b - 3c sin 20 , 16 

(a - 0)A . 
= {l + a 2 )A 2 [] ~ ]cSm{m > 

where the critical current j c is given by 



TrrA 



1 + a 2 



(1 + « 2 )A§- 



3c = 3A-. 



sinh(7iTA) : 



aKA 

\^W\ 



(17) 



(18) 
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FIG. 2: (color online) Velocities zo and <j> given by Eqs. (16) 
and (17) at j — l.lj c vs. time. (In arb. units) 



j'a is a critical current for the domain wall of the same 
width, but without the twist. These equations reduce to 
Eq. (13) for K = 0. We also note that Eq. (18) is correct 
only in the first order in K. 

The critical current j c is exponentially suppressed for 
twists TA > 1/-7T. For small twists j c w Ja(1 — 
7r 2 (rA) 2 /6). Note that j c in Eq. (18) diverges at a = 0, 
that is the DW does not spin for any current [17]. 

For j < j c the DW tilts by the angle sin(2(/>j) = j/j c 
and moves with a constant velocity zq = j/3/a, if = 0, 
the DW does not move at all. For j > j c the DW both 
spins and moves along the wire. 

Eqs. (17) and (16) can be solved analytically. The solu- 
tion gives both the velocity and angular velocity, which 
periodically depend on time [18] (see Fig. 2), with the 
period T and average angular velocity CI given by [11]: 



0=^ = 
T 



(a - 0)A 



(19) 



More experimentally relevant, however, is the average 
(drift) velocity of the DW Vd = (io)r- For any current 
it is given by [11] 



Yd 



for j < j c , 



I p . (a-gKl+ofA) m, 72 fnr „■ ^ a 

\aJ-\ a(l+ a *) VJ ~Jcj IOrj>J c . 



(20) 



The square of the deviation of the velocity from the drift 
velocity, Eq. (20), ((SV) 2 ) T is 



((6V) 2 ) T = 



0, 



q(1+q 2 ) J<- 



for j < jc 
for j > jc 



(21) 

Both Vd(j) and ((SV) 2 )t are shown in the inset of Fig. 1. 
For large currents, j j c , the drift velocity asymptot- 
ically approaches the velocity given by Eq. (13), while 
((8V) 2 )t approaches a constant. 

Summary. We have studied the effects of DMI on the 
magnetization statics and dynamics in a thin ferromag- 
netic wire. We have derived a simple criterion which de- 
termines whether the wire with the spiral magnetization 
state can sustain a DW configuration. Namely, in the 
wires with weak enough uniaxial anisotropy and/or ex- 
change constant compared to DMI constant (2JA < D 2 ) 



a DW cannot be formed. In the opposite case (2JA > 
D 2 ) we have found the spiral magnetization state with 
a DW in the wire. For = the wall moves along 
the wire only if the applied current is above j c given by 
Eq. (18). The variance of the velocity in this regime is 
given by ({5V) 2 ) T = V*(j/y/p=j* - 1). For jt 
the DW moves but does not rotate for currents below j c . 
Above j c the DW both moves and rotates [18]. Our re- 
sult, Eq. (18), shows that the critical value of current is 
suppressed by DMI. We also have derived the expression, 
Eq. (20), for the drift velocity Vd of the DW for all values 
of current. It shows that above the critical current j c the 
drift velocity can be enhanced by DMI. 

We believe that our findings can be experimentally ob- 
served, e.g., with the use of the scanning tunneling mi- 
croscopy which was employed to reveal the DW structure 
in ultrathin Fe nanowires [19, 20]. We note that in a real- 
istic experimental setting besides the "intrinsic" pinning 
there always going to be an extrinsic pinning due to a 
nonideal shape of the wire. It is, however, clear that in 
the near future the development of better nanofabrica- 
tion techniques will lead to the situation when one has 
to worry mostly about the "intrinsic" effect. 
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